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Abstract 



New extra series of conserved densities for the polytropic gas model and 
nonlinear elasticity equation are obtained without any references to the re- 
cursion operator or to the Lax operator formalism. Our method based on 
the utilization of the symmetry operators and allows us to obtain the densi- 
ties of arbitrary homogeneity dimensions. The non polynomial densities with 
. logarithmic behavior are presented as an example. The special attention is 

paid for the singular case (7 = 1) for which we found new non homogeneous 
■ solutions expressed in terms of the elementary functions. 
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CO 

Introduction 

The conservation laws are the most important object in the classical mechanics as 
well as in the field theory. There are many different methods of the constructions 
of these laws. The most popular, especially used in the soliton theory and in the 
hydrodynamics, utilize the so called recursion operator or Lax operator formalism 
PQ. On the other hand, it appeared that for the Nonlinear Schroedinger equation 
it is possible to find one series of conserved Hamiltonians using recursion operator 
only. However for the shallow water equation, which is the dispersionless limit of 
the Nonlinear Schroedinger equation, there is additional series of conserved densities 
which is impossible to obtain by recursion operator (see for example 0). 

In this paper we would like to show that it is possible to construct new extra 
series of conserved densities for the polytropic gas model and nonlinear elasticity 
equation [31 E] avoiding using recursion operator or Lax formalism. More precisely 
we generate many nonequivalent Hamiltonians of the given dimensions, using sym- 
metry operator. Our Hamiltonians are non polynomial expressions which contains 
logarithmic functions. Independently we consider the singular case, for the poly- 
tropic gas system for which (7 = 1) . For this system we constructed new series of 
non homogeneous solutions expressed in terms of the elementary functions. 
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The paper is organized as follows. In the first section we describe the basic 
properties of the polytropic gas system which we use in the next sections. In the 
second section we describe our symmetry approach where we utilize the shift, scaling 
and projective operators in order to generate the conserved densities. In the third 
section we present explicitly new series of the conserved densities for different, phys- 
ically interesting, models of the polytropic gas system with 7 = 2, 3, 4, 5, |, |, — 1. 
In the last section we adopt our formalism to the degenerated case (7 = 1) which is 
obtained from the contraction of the previous case. We show that in this case the 
conserved densities are connected with the Bessel equation. 



1 The hydro dynamical systems. 

The theory of the hydrodynamical type systems of the nonlinear equations [5] 

uj = X>j( u K i,j = l,2,..,N (1) 

i=i 

where u = (m, u 2 , ...un) and v* are some functions , integrable by the generalized 
hodograph method [H] is closely related to the over-determined systems of first order 
partial differential linear equations. The conservation laws are such that 

dh = dg_ 

dt dx 1 ' 

where h is density and g is flux. Then densities satisfy 

d k (d l h)=T\ k (d l h) + T k lk (d k h). i^k, (3) 

where 

T\ h = -^- i^k d k = d/dr k (4) 

and r fc (u) are the so called Riemann invariants in which the hydrodynamic type 
system (JTJ is rewritten in the diagonal form 

r\ = //(rK, (5) 

and no summation on the repeated indices. Thus, is a linear systems of first 
order partial differential equations with variable coefficients. The general solution 
of such system is determined up to N arbitrary functions of a single variable. 

There are many particular cases of the system Q for which a general solution is 
expressed in explicit and in compact form ( j7j ) . If we cannot to find a general solution 
of such system, then alternative way to solve the Cauchy or Goursat problems is to 
create the infinite number of particular solutions [Hj. 

It appeared that the conserved densities can be used in the construction of par- 
ticular solutions. Indeed. Let us consider the polytropic gas 

T)t = d x (w]), u t = d x — H , (6) 

L I 7 — 1 J 
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and nonlinear elasticity equations 

;„ 11... = rL 

W-2 



r v 1 ' 2 i 

T)y = U X , Uy = d X , (7) 



which are the commuting flows to each other. 

The solution of the first system obtained by hodograph method is 

t h ^ dh udh 

drj ' di] r] du 

while for the second is 

dh dh 
X dr]' ^ du 
where h is some solution of the Tricomi like equation 

Ku = ?7 3_7 /Vr ( 10 ) 

The previous equation is nothing but the equation on the conservation law densities 
for both systems. Formulas (jSJ) and realize the general or particular solution for 
both systems if h is the general or particular solution of the equation (fTUj) respec- 
tively. 

In jHj and jl] two infinite serieses of quasi linear conservation laws were con- 
structed for the polytropic gas and for the nonlinear elasticity equations respec- 
tively. Quasi linear means, that these conservation law densities are polynomials 
with respect to u, rj and rp where 7 is an arbitrary polytropic constant. 

However, as we see in the next section, these densities does not exhaust the set 
of all possible conserved densities. 

2 Symmetry operator approach for 7 ^ 1. 

Let us first consider the more general form of the polytropic gas system 

u 2 

T) t = d x (u V ), Ut = d x [-+ V f"( V )-f'( V )}, (11) 

Vy = u x, u y = d x f"(r}) (12) 

where / is an arbitrary function. Our generalized system constitute Hamiltonian 
equations with the following local structure 

„. a « Ut = d JJL, % ^f, ^f, (13) 

du Or} du or] 



where 



H = J[^u 2 ri + V f(r l )-2f{r ] )}dx, 



u 2 



H = H- + f(v)]dx. (14) 
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In order to obtain the conserved densities we try to eliminate the differentials dq 
and dp from the corresponding conservation laws 

d t h = d x p, d y h = d x q (15) 

By direct calculation we have 

dq = h^du + f'"(r])h u dr], dp = [uh u + r/h^du + [i]f"\ri)h u + uh v ]dr] (16) 

The compatibility conditions ( (p u ) v = (p v )u, (<?«)»? = (qr)u ) lead to the Tricomi - 
like equation 

Kv = f m {v)Ku, (17) 

Since, this equation is compatible with the shift symmetry operator (see [TD1 ) 

5 = d/du, (18) 

one can search solutions in the form 

K = Xh, (19) 

where A is an arbitrary parameter. This is the eigenvalue problem for the shift 
symmetry operator. Then Tricomi-like equation is transformed to the linear ordinary 
differential equation 

h m = X 2 f'"(v)h- (20) 

However this equation cannot be solved explicitly for an arbitrary /. 

We propose, from that reason, quite different approach. Let us observe that the 
shift symmetry operator transforms one solution of the Tricomi - like equation onto 
another one 

d u h n+1 = h n . (21) 

It means that all conservation law densities hk and corresponding fluxes pk,qk can 
be written down in the quadratures recursively 

dh k+ i = h k du + q k dr], dq k+1 = q k du + f"'(rj)h k drj, (22) 
dpk+i = [uh k + r)Qk\du + [i]f"'(v) h k + uq k ]drj. (23) 

Interestingly all known examples of the conserved densities could be obtained in this 
way. 

However we demonstrate different possibilities which give us new serieses of the 
conserved densities. Our main idea is based on the following observation. 

The Tricomi like equation (fTU|) is compatible with three local symmetry opera- 
tors: 



1. Shift operator 

2. Scaling operator 



s =l < 24 > 



d 2d 

R = U 7T + T^TT' ( 25 ) 

ou 7 — 1 or) 



3. Projective operator 

S = + ( 7 - l)- 1 ^ 1 ]^ + + l^u. (26) 

These operators, with the identity operator, constitute closed Lie Algebra, with the 
following commutation relation 

[<5,S] = ^i?+^, [S,R}=5, [R,S] = S (27) 

They act on homogeneous conservation law densities as follows 

Shk+i = hk Rhk = Ckhk Shk = hk+i, (28) 

where Ck are degree of homogeneity. Thus, by combination of these symmetry op- 
erators one can describe all quasi linear conservation law densities. 

Interestingly if we rewrite the Tricomi -like equation using the Riemann invari- 
ants we obtain the famous Euler - Darboux - Poisson equation 



dh e v dh dh 



dridr 2 r\ — r 2 ^dr\ dr 2 

where 



(29) 



£ = W^iY ri = u + ^h) v ^' r2 = u -^ L Y) 71 ^- (30) 

Now the symmetry operators are 

d d i ^ 2 ^ 

dr 1 dr 2 ' ^ dr 1 ~^ dr 2 ' ^ ^ 

S = (rT^r + ^^ + ^ + r 2 ]. (32) 

The "zero-solutions" of the shift operator 5 

5h = (33) 

can be found immediately 

= h^= V , (34) 
while for the projective operator one can obtain 

= u, h<i ] = urj, (35) 

etc. 

The "zero-solutions" of projective operator are easy to obtain also from Riemann 
invariant (r 1 , r 2 ) form 

h = (rV 2 )- £ . (36) 

If we shift the Riemann invariants, in this formula, on arbitrary parameter A we 
obtain the generating function on the conservation law densities. When A — > oo we 
can obtain known densities. 
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In order to obtain new densities we demonstrate quite different possibilities . Let 
us consider the following recursive chain 

Rhf = c k hf, RhP = c k [h^+hf\ Rh%> = h£>] (37) 

etc., where h% are the quasi linear conserved densities. 

Now h k 1 , % = 1,2, ... are new conserved densities which satisfy the Tricomi-like 
equation (jlUj) . Notice that the second formula in (|2*Hj) is not in contradiction to 
previous chain because formula (J2~%j) is valid for homogeneous conservation laws 
only. 

Now let us consider the shallow water equation, e.g. the polytropic gas system 
with 7 = 2. Let us choose the first homogeneous solutions of 



hf ] =U, h^= V , hf=UT) =U 2 7 ] + 7 1 2 (38) 

etc. As the result we obtained that 

= 2\Ju 2 - Ar][hi (■\Ju 2 - 4i] + uj - In 2) + (u - \ju 2 - 4rfj In 77 

, (1) u2 

h 2 = — + r] In T] 

h*p = - (u 3 + 6ui] In 7] 



h 



(i) 



\(u A + I2u 2 rj\nrj + Urflnrj- 18r? 2 ) (39) 



are conserved density also. In the next let us apply this procedure to the Ivp and 



h {2) = — f-36r]ln 2 (^/z + u) + 18ln(y/z + u)(2T]uln(Ar])-2r] + uVz) 
18 ^ 

+9uyG(l - ln(4r?)) - 36r]lnr/ln2 + (9m 2 + 18t?) lnr/ - 18m 2 + 32rj). 



h-p = — ( — 36T]u\ri 2 (y/z + u) + ln(\/z + m) (36t7m ln(4r^) — 6i]u + 

&\fz{u 2 + 8r])) — 36r]u ln(^) In 2 + 3 lnr^u 3 + 6r]u — \fziv 2 + 8r/)) + 
3-^(1 -21n2)(w 2 + 877) - I6r]u-8u 3 ), (40) 

where z = u 2 — 4i]. 

We obtained rather complicated formulas which contain the logarithm in the 
second power. Now we can continue this procedure obtaining complicated formulas 
which contain the logarithm in higher powers. On the other side we can combine our 
method with this which generate quasi linear densities obtaining following diagram 
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R 



h, 1 



R 



h,° 



k+1 



R 



h! 



k+1 



R 



h 



k+1 



This diagram obviously is closed. In order to see it let us start consideration from 
some h k . Applying at first projective and at next step scaling transformations we 
obtain 

Rhft? = c k+ i(ht\ 1] + ShP) (41) 

Finally applying these operators in the reverse order we obtain 

SRh ( k n+1) = c k (Sh[ n+1) + Sh[ n) ) (42) 

Using the commutation relation (|27j) we see that these formulas coincide because 
c k+ i = c fc + 1. 

3 Non polynomial densities for an arbitrary 7 7^ 1 
in polytropic gas system 

For an arbitrary 7 our generating equations become 

H4 n) = c fc (4 n) +4 n_1) )- («) 

(n) 

Now we have to solve the Tricomi like equation on h k which is however rather hard 
to do. We restrict, from that reason, to two cases where hf 1 = u and — V- 
For the first case we have C\ — 1. By direct verification one can check that 

= uh(r) + u\nu (44) 

where r = satisfy (|43p. The Tricomi like equation reduces in this case to 

d 2 h , ^ \dh 



>j 2 

)2l 



r 2 ' 



(4 " (7 - I) 2 ) ^ + (6 - 7(7 - l)r) ^ + 1 = 0. (45) 



For the second case we have C2 = Similarly to the first case, by direct 

verification, one can check that 

= f]h(r) + i] In r) 

(46) 
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satisfy ()43j). Substituting to the Tricomi like equation (fTUj) we obtain the fol- 
lowing equation on the function h(r) 

r) 2 h f)h 

r(4 - r( 7 - 1))^ - (2 - (7 - - 2)t)^ + 1 = (47) 

In both cases it is possible to obtain closed formulas for particular values of 7 
parameter. However we can simplify our consideration using different coordinates. 
For example choosing the coordinates r and p as [3] 

r] = rp, u = (r 7_1 + p 7_1 ) , (48) 



7-1 

we can rewrite the scaling symmetry operator and Tricomi like equation as 

R — (r— H —1 

7 — 1 \ dr ^dp' 

P'-^ p = ^ (49) 

Taking into an account hf^ = r we obtained the following solutions on the hi 
densities 

7 = 2, /i^ = — (r — p) ln(r — p) + p\n(p) 

7 = 3, = (r + p) ln(r + p) — (r — p) ln(r — p) 

7 = 4, = — (r — p) ln(r — p) + (r + 2p) ln(p 2 + pr + r 2 ) 



1 p + 2r x 

+V 3parctan( — = — 
V3p ' 



7 = 5, /ij = (p + r) ln(p + r) + (r — p) ln(r — p) + r ln(p 2 + r 2 ) 

+2parctan(^) (50) 

Finally let us present the solution for 7 = |, | and 7 = — 1 which are interesting 
from the physical point because these describe dynamics of one-atomic gas and 
two-atomic gas [llj and two dimensional nonlinear Born-Infeld electrodynamics [T2"] 
respectively. 

5 qi 3 1 t , 2 iox 

7 = -, — -p arctanhrs + - In (3r» — 3rs + r — 1) 

_ 2 

— rmr — 3rr 3 

T 5 7* 

7 = -, h\ = -p arctanhrs -\ — In ( — 10rs — 5rs + 5rs + r= + r 2 — 1) 

5 4 2 

1 5 _2 _4 

— r In r — -7T 5 — 5rr 5 + r In r 
7 = -l, h? = Ur\n(^- 2 )+p\nC-^)) (51) 
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4 Symmetry operator approach for 7=1 and 
corresponding non polynomial densities. 

For this case the polytropic gas and the nonlinear elasticity systems are 

vU 2 1 

ilt = d x (ur]), Ut = d x [— + In 77J . (52) 
rj y = u x , u y = d x ( - i), (53) 



respectively. For this case the Euler - Darboux - Poisson equation (jzSj) degener- 
ates because e — > oo. Similarly the scaling operator R (see (f23|) ) also becomes 
degenerated. 

In order to solve this problem we notice that the analogue of the symmetry 
algebra ()26|) can be obtained considering the contraction with respect to 7 = 1. If 



we rescale R — ► ^-R and compute the limit when 7 — > 1 we obtain 

5 = d u , R = r]d v , S = \nrid u + urjd v — u/2, (54) 



Interestingly now, in the contraction limit, the homogeneity properties ()28|) does 
not hold. It means that the Tricomi-like equation 

h U u = V (55) 

has no any homogeneous solutions. 

However the quasi- homogeneous solutions could be obtained ( 9j). By quasi- 
homogeneous we understand the homogeneous with respect to the one variable u or 
In 77. 

For our further applications we describe another possibility. We constructed 
the generating function for the conserved densities. To end this we consider the 
eigenvalue problem for each symmetry operator. 

1. For the scaling operator R (see (154)0 the eigenvalue problem 

hz = Xh, (56) 

where = r]d v , reduces the Tricomi-like equation to the linear ordinary differential 
equation of the second order 

h uu = A(A - l)h. (57) 
These equations can be easily integrated and we obtained 



/i = exp[yA(A-]> + A£]. (58) 

2. For the shift operator 5 (see 1)540) the eigenvalue problem 

h u = Xh (59) 

reduces the Tricomi-like equation to the linear ordinary differential equation of the 
second order 

fc K -ht = X 2 h, (60) 
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These equation are easily also to integrate and the solutions have the same form as 
()58j) in which we should replace A 2 — > A(A — 1). 

If one expands function h near A — > or A — > 1, then one can obtain two 
well-known infinite serieses of conservation law densities (0). 

3. Projective operator S (see (JSH))- If one substitutes 

h = qexp[£/2], u = 2s cosh 9, £ = 2ssinh6 l , (61) 
where s and 9 are new functions, then the eigenvalue problem looks like 

qe = Xq (62) 

Now the generating function is 

q = ip(s) exp[\9], (63) 
where the function if) satisfies the Bessel equation 

if)" + V + [1 - K]*i> = 0, (64) 
s s z 

Let us present the most simplest example of the conserved densities obtained 
in this way, where it is possible to describe the Bessel function by the elementary 
functions for which A = | 

Using the recursion relations, known in the theory of Bessel functions, it is possible to 
generate infinite sets of conserved densities written in terms of elementary functions 
without any references to the recursion operator, which appear in the Hamiltonian 
approach to the polytropic gas systems. Finally we can use the shift or scaling 
symmetry operators and generate some conserved densities as we did in the previous 
sections. 



5 Conclusion 

In this paper we constructed new extra series of conserved densities for the polytropic 
gas model and nonlinear elasticity equation avoiding using recursion operator or Lax 
formalism. Our Hamiltonians appeared as the non polynomial expressions which 
contain logarithmic functions. If we continue our procedure to these logarithmic 
densities in the next step we obtain expressions with the logarithmic function in 
an arbitrary power. We considered the singular case of the polytropic gas system 
also, for which we found the non homogeneous solutions expressed in terms of Bessel 
functions. If we apply the point transformation to the symmetry operators then they 
can change the role. It means that the scaling and shift operators are equivalent 
to each other, what it easy to see in the Riemann invariants ()31j) or (J59I56J) . We 
presented this method for two hydrodynamical systems only, but this method can 
be adopted to more complicated equations as well. 
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